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Geometric phases arise in a number of physical situations and often lead to systematic shifts in 
frequencies or phases measured in precision experiments. We describe, by working through some 
simple examples, a method to calculate geometric phases that relies only on standard quantum 
mechanical calculations of energy level shifts. This connection between geometric phases and off- 
resonant energy shifts simplifies calculation of the effect in certain situations. More generally, it also 
provides a way to understand and calculate geometric phases without recourse to topology. 
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I. INTRODUCTION 

A geometric phase, often also referred to as an adiabatic phase or Berry's phase, is a real, physical phase shift that 
can lead to measurable effects in experiments ([IJEIIB] are some examples). In the classic version of this effect, the 
state of a particle with a magnetic moment is modified under the influence of a magnetic field which undergoes a slow 
change in its direction. The motion of the field leads to a phase shift accumulated between the quantum states of 
the particle, in addition to the dynamical phase due to the Larmor precession of the magnetic moment around the 
i— | magnetic field. This extra phase is termed a geometric phase because it has a simple interpretation in terms of the 
geometry traced out by the system's Hamiltonian as it evolves in its parameter space. We point the interested reader 
Q-i to [1] and references within for an excellent account of the geometric approach to calculating geometric phases, and 
its various connections with the topology of the parameter space. 

In practice however, the formulation of the geometric phase in terms of differential geometry can make its calculation 
somewhat cumbersome. For example, calculating geometric phases can be non-trivial when the path traced out in 
the parameter space is not a closed loop [5]. In addition the geometric formulation is strictly correct only in the 
i— i limit of adiabatic (slow) evolution, and corrections due to finite evolution speed can be important in some situations 
[6 . Finally, the parameter space of the Hamiltonian can also become quite involved when the system is an atom 
or molecule with internal structure subject to various evolving fields (e.g. magnetic fields, electric fields, and their 
gradients). Far from being an abstract intellectual pursuit, our interest in calculating geometric phases for this sort 
of complicated system stems from precision atomic physics measurements, such as the search for permanent electric 
i-H dipole moments [7], where these phases can give rise to important systematic errors [TJ [5J [5] and/or quantum state 
decoherence [Ttjl ITT] . 

^ i Energy level shifts from time-varying perturbations, such as Bloch-Siegert shifts and AC Stark shifts, are mostly 

treated in the literature as separate phenomena from geometric phases. However, the idea that both of these essentially 
involve the same physics is implicit in some recent work [3i| 8 , 9, 12 and originally derives from an analysis by Ramsey 
|13j . While most of these treatments analyze the simple case of a spin- 1/2 particle in a time-dependent magnetic 
field, the link between geometric phases and off-resonant energy level shifts is in fact quite general. In this paper, 
we explore this connection by looking at some instructive examples, showing with simple perturbative calculations 
how energy level shifts lead to geometric phases in each case. We find that there are a number of advantages to this 
reformulation of the problem of calculating geometric phases. The energy shift approach to geometric phases can be 
applied to situations where the geometric formalism is opaque, and our calculations of energy shifts using perturbation 
theory only rely on standard undergraduate quantum mechanics. We feel that it is pedagogically useful to highlight 
the equivalence between geometric phases and the more familiar energy level shifts, thereby showing that there is 
nothing arcane about geometric phases. 

In the following section we introduce our approach with the classic example of a spin-1/2 system subject to a 
magnetic field whose direction changes in time. Then in Section [III| we work out modifications to the usual geometric 
phase when there are multiple Fourier components in the time evolution of the magnetic field. Sections |IV| and [V] 
show how the geometric phase can arise in a system acted upon by electric fields or a combination of electric and 
magnetic fields. These sections also indicate how geometric phases can be analyzed in a system with a non-trivial 
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FIG. 1: An example of a trajectory traced by the magnetic field, showing how it can be decomposed into its components along 
a fixed set of axes: B = B z z + B±(t). The geometric phase <j> g picked up between the sublevels of a spin-1/2 system interacting 
with this field is equal to the solid angle fi s enclosed by the trajectory. 

level structure, such as an atom or molecule. 

II. SPIN-1/2 SYSTEM IN A MAGNETIC FIELD 

A spin S system in a magnetic field is a well-studied example of the geometric phase [TJ 13 H]- When the tip 
of the magnetic field vector slowly traces out a closed loop in space, the sublevels \ms) of the system pick up an 
extra "geometric" phase <j> g (ms) = rns x fl s (in addition to the usual dynamical phase from Larmor precession) 
[U Q31 Q1)J US]. Here f2 s is the solid angle enclosed by the loop. In all the subsequent sections, without any loss of 
generality we focus on the observable quantity defined by the phase difference between the \ms = ±<S) sublevels: 
A(j) g = <f>{ms = +S) — 4>(ms = —S). In this section we will examine a spin S — 1/2 system in a magnetic field 
and show how the geometric phase shift can be traced back to an off-resonant energy level shift, in this case an AC 
Zeeman (or Ramsey-Bloch-Siegert) shift. 

The usual approach to calculating the geometric phase tracks the instantaneous eigenstates of the system under 
the assumption of adiabatic evolution j4j HH [15] [16] . In contrast, we shall use a fixed coordinate system and resolve 
the changing magnetic field into a static longitudinal (z-directed) component and a dynamic transverse component 
(in the xy plane), as shown in Figure [I] Let the evolution of the transverse component's direction be composed of a 
single angular frequency luj_. For a counter-clockwise rotation, the magnetic field written in terms of its components 
in the fixed coordinate system is 

B = B z z + B ± (t) (1) 
where the rotating transverse component of the magnetic field is 

B±(t) — B±(x cosuj±t + y sinwjj). (2) 

B 2 

For this discussion we consider small values of Bj_/B z , such that the solid angle enclosed by this loop is fl s w it -g^. 
The Hamiltonian of the particle in this magnetic field is 



Hint = -jS ■ B 



(3) 
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FIG. 2: Energy levels of a spin S — 1/2 system interacting with a magnetic field whose direction evolves over time. The level 
splitting in the longitudinal field is hu>o — yB z . The dashed lines indicate virtual levels, offset from the real levels by the finite 
rotation frequency of the transverse field B±. S± is the spin raising (lowering) operator which couples the \ms = ±1/2} state 
to the [virtual] \ms — Tl/2) state in the presence of B±. 



= -jS z B z 



— iuj±_t 



s 4 



(4) 



where 7 is the gyromagnetic ratio and S± = 
presence of only the longitudinal field B z z, 



- S x ± iS v are the spin raising and lowering operators respectively. In the 
the eigenstates \ms = ±1/2) have energies E ms = (ms\ — ^SzB^ms) — 



—^B z m s and the Larmor precession frequency is luq = (E_ 1 / 2 — E +1 / 2 )/Ti = ^B z /fi. We will consider the effect of the 
time- varying transverse field as a perturbation on the S z eigenstates. To lowest order, the energy shifts A_E ±1 / 2 of the 
I ± 1/2) states are given by second order time-dependent perturbation theory. Figure [2] shows the energy levels and 
the operators connecting them. We define the transverse matrix elements of the spin operator to simplify notation 
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|<+1/2|S + |-1/2)| 2 = |(-1/2|S_| + 1/2)| 5 



1. 



(5) 



Taking note of the usual selection rules on the matrix elements of the S± operators, we get 
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l 2 s 2 ± Bi 
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B z 



(7) 



Here we have retained terms up to lowest order in u±, the field's evolution frequency. Note that the condition for 
applicability of this approximation viz. luj_ <C uiq is the same as the adiabatic criterion invoked in the standard 
approach to describing the geometric phase |14j . 

The extra energy difference AE between | ± 1/2) due to the transverse rotating component of the field is therefore 



AE = AE + x/2 — AE_i/ 2 

r B\, . is? . 



AE, 



QZ 



AE, 
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(8) 
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where AEqz (AE g ) corresponds to the term in curly (square) brackets. Note that the term AEqz is nonzero 
even when the field's direction evolves infinitesimally slowly. This term is nothing other than the correction to the 
Zeeman splitting because the total magnetic field becomes larger on application of the transverse field component: 

here B = y'BJ + £> 2 _ « B z + The term AE g is responsible for the geometric phase. This term, which vanishes 
when luj_ — *• 0, nevertheless adds a relative phase between | ± 1/2) even in this limit. The relative phase picked up 
over a time interval T is 

^ g = -AE g T/?i = ^u ± T. (9) 

Over one full cycle of the field's evolution, T = 2tt/lu± and we get the standard geometric phase result A(f> g = 
B 2 

TT-gt — The restriction to small solid angles here is equivalent to truncating the perturbation scries at second order 

z 

in B± . Results in the case of larger solid angles can be analytically written down using higher orders of perturbation 
theory, or alternatively using a dressed-state formalism that is valid for arbitrarily large solid angle |12l I17j . However 
the essential point remains that energy level shifts can be calculated algebraically (by numerically diagonalizing the 
Hamiltonian if necessary), without having to track the eigenstates over a geometric path in parameter space. 

This formulation can be easily extended to situations that are more complicated than the simple loop. Moreover, 
by keeping terms of higher order in u>±, corrections due to deviations from the adiabatic limit can be calculated. 
The energy shift formalism also eliminates the need to track the instantaneous quantization axis and basis states of 
the system throughout its evolution. The use of a fixed basis set makes contact with the standard tools of quantum 
mechanics as taught in most undergraduate courses. Having illustrated the basic idea, we now apply this formalism 
to a case where the evolution of the field's direction is more complicated. 



III. SPIN-1/2 SYSTEM WITH MULTIPLE EVOLUTION FREQUENCIES 
A. Geometric phase calculation using energy shifts 



Let us calculate the geometric phase in a situation where there are two Fourier components to the evolution of 
the magnetic field experienced by a spin S — 1/2 system. This demonstrates the approach to be used in a general 
situation, e.g. in a precision experiment where the magnetic field experienced by an atom along its trajectory can be 
quite complicated and contain a number of Fourier components. A magnetic field with two evolution frequencies u>\, 
lu 2 is 



B = B z z + B±(t) 
B±(t) = B\{x cosuit + y sinujit) + Bi{x cosuj 2 t + y sinc^i) 

Using the same analysis and notation as the previous section, the energy shift of the ground state is 



AE +1/2 = 



1 



huji — huj hu 2 — hujo 



l 



Uui 2 — Uojq 



(10) 



(11) 



where we have retained the terms in second order perturbation theory at the beat frequencies between the two Fourier 
components of the field. The energy shift of the excited state is 
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—hui 2 + huiQ 



(12) 



Define the difference between the two evolution frequencies to be Aw = u)% — u)\. The extra energy difference 
between the excited and ground states due to the transverse magnetic field is AE = AE +1 / 2 — AE_ 1 / 2 . Up to first 
order in the evolution frequencies it is given by 



AE = 



'7 



2B Z 
B\ 



B 2 
2B Z 



6162 a A 
7—^ — cos Aujt > 



B z 



B1B2 
26? 



h(u 2 + uj\) cos Aut 



AE QZ + AE g , 



(13) 
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where AEqz (AE g ) corresponds to the term in curly (square) brackets. 

All the oscillating terms have been retained above. The phase difference arising from the geometric term is 



A<j) g = / -AE g dt/h = 



1 



B\ 



B 2 

UJ\T + ~rp:i0 2 T 



B\B 2 (ui + lu 2 ) 



Bl 



Ao 



(14) 



Note that here the geometric phase is not just the sum of phases due to the fields B\ and B 2 , but also contains an 
oscillating cross term that is proportional to sin AloT. This term depends on the phase relationship between the two 
transverse field components, and cancels (or not) depending on the difference frequency Alu and the evolution time 
T. This is a physically measurable effect that arises from time-dependent energy shifts that are usually assumed to 
average away. In the following section, we verify this result by calculating the effect of two evolution frequencies using 
the usual geometric method. 



B. Geometric approach 



For a spin S =1/2 system, the geometric phase difference is Acf> g — £l s . We will evaluate the solid angle f2 s enclosed 
by the tip of the magnetic field vector defined by Equation (10 1. The solid angle enclosed by a curve on a sphere is 
given by 



n s = J (1 - cos 9) dcj> 

where 0, <j> are the usual spherical polar angles. We use the following definitions 

B = yjBl + Bl 
B z = Bcos9 



(15) 



(16) 



and find that the solid angle is 



n s = 



B-B z 

B 
1*1 
2 Bl 



A 

Bl 



(17) 



where A = J \B\ dcf) is the area enclosed by the transverse field B±(t) in the ccj/-plane. For the purpose of comparison 
with second order perturbation theory the approximation above, where we retain terms up to second order in B±, is 
sufficiently accurate. 



We can express the area A in terms of B±(t), and evaluate it for the curve defined in Equation ( 10 1: 

1 

,u 2 



A 



1 



-Bl dcj) : 



B±(t) x 



dB±_ 

~dT 



■ z dt 



1 



dt 



{B\ cosujit + B 2 cos uj?£)(w\B\ cos u)\t + uj 2 B 2 cos u 2 t) 



(Bi sin wit + B 2 sin u>2t){u>iBi smujit + lo 2 B 2 sina^) 



1 
1 



dt 



B\liJi + B\w 2 + BiB 2 (loi + w 2 ) cos Awi 



BluoiT + B\w 2 T + BiB 2 ^— A + ^ sin AuT 

Auj 



(18) 



Using Equation (17), we find that the geometric phase difference is 

A 



A(j) g = Q s = 



B 2 



^\B{ rp.Bi rp , BiB 2 (wi + w 2 ) . 

^T + -u 2 T + — ^— sin AuiT 



(19) 
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exactly as in the calculation using perturbation theory in Equation (14 1. We note that Ref . jlj discusses a particular 



example of this situation, where u>\ — 2cl>2 and hence the tip of the magnetic field vector executes a figure-8 motion. 

IV. SPIN-1 SYSTEM IN AN ELECTRIC FIELD 

Using the same formalism as above, the geometric phase can be calculated for a system with more complicated 
levels such as an atom or a molecule. In this section we will consider the geometric phase for a spin J = 1 system 
in an electric field whose direction changes with time. As the electric field vector traces out a loop, the system picks 
up a geometric phase A(j) g = 2fi s between the \mj — ±1) states. We refer the reader to the calculation for this case 
using the geometric formalism in 1 16] . and calculate the same here using energy shifts. 

Let the evolving electric field be written as 

£ = £ z z + £ ± (t) (20) 
£±(t) — £±(x cosujj_t + y siucjjj) (21) 
= ^=(-f +1 e-^* +f_ 1 e+^ t ) (22) 

where r±i = T ^ X ^P are spherical basis vectors. 

The Hamiltonian of the system in the electric field is 

Hint = —D ■ £ 

= -D z £ z ^ +1 ± . (23) 

Here we have defined the electric dipolc moment operator D. This operator only couples states of opposite parity, 
whereas the \mj — 0, ±1) sublevels of a J = 1 level all have the same parity. To calculate the effect of electric fields, it 
is essential to enlarge the system and include an opposite parity state in addition to the spin-1 sublevels. For simplicity, 
we consider here the 4-state system consisting of the 3 sublevels of a J v = 1~ level: | J = 1, to,/ = 0), | J = 1, mj — ±1), 
and in addition a J x = + state: | J = 0, mj = 0). Here J denotes the angular momentum and it the parity of the 
state. We refer to states by their \J,mj) labels from now on. Figure 3]ja) shows these states. We are interested in 
the phase that is picked up between the |1,±1) states due to the evolving electric field. Choose the zero of energy 
halfway between |1, 0) and |0, 0), and let the zero-field separation between them be 2B. The perturbative calculation 
requires matrix elements of D, which we write as 



d z z = (1,0|£>|0,0) 
-d ±1 f Tl = (1,±1|£|0,0>. (24) 

Here we have used the Wigner-Eckart theorem to define only the nonzero matrix elements of D; this can also be used 
to show that all three of the nonzero matrix elements have a common value: do = d z = d±\. 

First, we consider the effect of £ z . This part of £ only mixes the two states |1,0) and |0,0). Under the interaction 
with this part of the field, the eigenstates are 

| CO) = 1 1,0) cos e/2 + |0,0) sine/2 

|070) = |0,0)cos£/2- |l,0)sin£/2 (25) 
where we denote the field-mixed eigenstates with tildes, and the mixing angle is 

d a £ z 



tan£ = 
sin 2 £/2 = 



B 



y/ff 2 + riggf - B 

2 V / B 2 ~Tdp2 



cos^/2 = ^] ^+ B . (26) 
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FIG. 3: Energy levels of a spin J = 1 system in an electric field. States are labelled by their angular momentum quantum 
numbers J,rnj. Dashed arrows and their labels indicate the states being coupled by the respective components of the dipole 
moment operator. Under the influence of a longitudinal electric field £ z , the |0,0) and |f,0) states shown in (a) are perturbed 
into 1 0,0} and |f,0) as shown in (b). The rotating transverse electric field £±{t) induces virtual energy levels (dashed lines) 
huj± above and below the static perturbed states. £±(t) couples these states with the |1,±1) states of interest, through the 
dipole moment operators D±i. 



The |1,0), 1 0, 0} states have energies ±^ B 2 + d%£ 2 respectively. The |1,±1) states remain at their zero-field location 

= -Eji,-!) = +B, as shown in Figure 3[b). The splitting Ai between the |1,±1) and |1,0) levels is usually 
referred to in the literature as a tensor Stark shift |16j . 

Next we consider the interaction with the rotating transverse field and calculate the effect to second order in 
perturbation theory. As indicated in Figure [ijb) the transverse field only couples |1,±1) with |1,0) or |0,0). Define 
the following energy denominators for ease of notation 



A = E^ ±1) -E m =B+^B^ + dl£^ 
A! = 25| 1)±1> - Epfi = B — ^B* + d 2 £ 2 . 



(27) 



Note that as defined, Ao (Ai) is positive (negative). The energy shifts AE±i for the |1,±1) states due to the 
transverse electric field £± are given by 



AE +1 = 
A£_i = 



dl£l cos 2 f/2 dl£lsm 2 £,/2 



Ao + huj± Ai + ~hoj± 
dl£l cos 2 Z/2 d^£l sin 2 £,/ 2 



Aq — Hu!j_ Ai — hbJj_ 

Where the |1, ±1) states were degenerate before, they are now split by the energy AE given by 

AE = AE +1 - A£_! 

AZ&rnJt/O d 2 £ism 2 Z/2^ <n/ 2N 



(28) 



A 2 



A? 



(29) 



We have again retained only the terms to least order in luj_ to illustrate the adiabatic part of the phase. Substituting 
the values of £, A , Ai we get 



AE 



d 2 F 2 



y/B 2 + d 2 £l + B 1 

2^WTd 2 £ 2 
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^B 2 + d 2 £ 2 z - B 

£1 



b 2 + d\ei - b 



huj± 



£? 



TlU! I 



AE n 



(30) 



(31) 



Compared to the case with a magnetic field, in a pure electric field the only energy difference between |1,±1) is 
the geometric contribution AE g . The relative phase between |1,±1) after a time duration corresponding to a single 
complete cycle of evolution, T = 2ir/ujj_, is 



A(b g = -AE g T/h = %uj x T 
£ 2 



(32) 



and the phase shift again has a simple geometric interpretation in terms of the solid angle enclosed by the evolving 
electric field. 



V. COMBINED EFFECT OF ELECTRIC AND MAGNETIC FIELDS 



The formalism developed in the previous section can be easily extended to a more complicated case, where the 
system interacts with both electric and magnetic fields. The standard geometric approach in this case requires us to 
calculate all instantaneous eigenvectors of the system, and their gradients with respect to the 6-dimensional parameter 
space that describes the Hamiltonian and its evolution (3 degrees of freedom each for the electric and magnetic fields). 
Energy shifts on the other hand can be calculated in a straightforward manner, as we shall see in the following. 



A. Rotating electric field, static magnetic field 

We will first examine the same spin-1 system as in the previous section and understand what happens to it when 
a static magnetic field B — B z z is imposed along with the revolving electric field. The interaction Hamiltonian of the 
system is now 



Hi; 



-D-S-p-B 

-D z £ z - [i z B z 



72 



(33) 



with jl the magnetic moment of the J = 1 state. 

To illustrate the basic ideas we again use a perturbative approach, to study the influence of the transverse time- 
dependent components of £ on the eigenstates in the presence of only £ z and B z . The eigenstates of the system are 
shown in Figure |4j We will carry over the notation and definitions from the previous section. 

The essential effect of the B z field is in the energy denominators that appear in the AC Stark shift. The energy 
shifts of the |1, ±1) states (in addition to their Zeeman shifts in the magnetic field) are given to lowest order in £± by 



AE, 



l r 

2 
1 



rfgg 2 cos 2 g/2 
A + \x z B z + ?iu)j_ 

d 2 £l cos 2 £/2 
A - [i z B z - Tiuj± 



d 2 £\ sin 2 C/2 
Ai + [i z B z + Tiuj± 
d 2 £\ sin 2 f/2 



Ai - UzB, 



(34) 



To make a connection with the adiabatic limit, we expand the energy shift AE between |1, ±1) up to first order in 
Tilu±, and write 



AE = AE +1 - AE-x kC + V ?uj±. 



(35) 



We interpret C as the analog of a (static) quadratic Zccman/Stark shift due to the inclusion of £±, and AE g = T> fruj± 
as the energy shift that leads to a geometric phase difference A(j) g = —Vlu±T. Tablc|T]shows the values of the geometric 
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FIG. 4: Energy levels of a spin J = 1 system in a combination of electric and magnetic fields. The state labels and symbols are 
the same as in Figure [3] In (b), the fi gure as shown corresponds to a case where the Zeeman shift is somewhat smaller than 
the tensor Stark shift: /j, z B z < Ai. 



phase A(f) g in some limiting cases, corresponding to regimes where one of the various energy scales in the problem is 
dominant. The relevant energy scales are the Zeeman shift fi z B z , the Stark matrix element do£ Zl the Stark shift of 
the 1 1,0) state |Ai|, and the zero-field energy splitting 2B between J 71 = + and J v = 1~ levels. 

Note that in general the geometric phase cannot be related to an obvious solid angle, due to the complicated 
parameter space of the Hamiltonian. While the evaluation of the actual solid angle by geometric means can be 
rather arduous, the geometric phase in every case can be simply calculated from the energy shift between |1,±1). 
This energy shift can be computed perturbatively, or by numerically diagonalizing the Hamiltonian if necessary. An 
important feature to note in Table [i] is that the zero-field energy splitting between J v = + and J" = 1~ levels 
appears irreducibly in some of the limiting expressions. This shows that the inclusion of the opposite parity {e.g. 
J n = + ) state is necessary to properly describe geometric phases in a combination of electric and magnetic fields. 



Case Limiting condition Geometric phase difference, A(f} g 



I ^B a <|Ai|,do£* 

II | Ai| < fi z B z < d £ z 

III do£z « VzBz < IB 

IV d S z , 2B < y, z B z 




TABLE I: Analytical expressions in some limiting conditions, for the geometric phase in a static magnetic field B z and electric 
field £ z , along with a rotating transverse electric field component £±. In each case, the geometric phase between |1, ±1} up to 
first order in ui± is shown. 
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B. Rotating magnetic field, static electric field 



As another illustration of the energy shift formalism in a case with combined electric and magnetic fields, consider 
the case where the system experiences a static i-directed electric field, and a revolving magnetic field B = B z z + 
B±(x cosuj±t + y smuj±t). The interaction Hamiltonian is 



Hint = ~D ■£ - Jx-B 
= -D z £ z - n z B z 



H_ x B^e~ iw ^ - n +1 B±e iu; ± t 

71 ' 



(36) 



where the spherical components of the magnetic moment fl are defined as follows. 



[i z Z 



(M|m|M> 

(i,±i|m|i,o). 



(37) 



Once again, the Wigner-Eckart theorem has been used to define the non-zero matrix elements (which have a common 
value): /i = fx z = fi±i. 
The energy shifts are 



AE +1 



2 



^icos^/2 



1 



Ai + n z B z + Uu!j_ 

nlB 2 ± cos 2 £/2 
Ai - \i z B z - hu± 



-A + \i x B z + hwx 
' /i^lsin 2 £/2 
-A - [i z B z - hujj_ 

As in the previous section, we expand the energy shift AE between |1,±1) up to first order in foui±_ 

AE = AE +1 - AE_ X kC + V huj x 



(38) 



(39) 



and focus our attention on the relative energy shift AE g = V Tiuj± that leads to a geometric phase A<f> g . The relevant 
energy scales in the problem are the same as in the previous section. We consider cases where Ai <C 2B i.e. the 
system is weakly electrically polarized by the £ z -field. Table [IT| lists the geometric phase in some limiting cases where 
the analytic expression becomes simple. 

As expected, when the Zccman splitting is sufficiently large (Cases I and II) the geometric phase reduces to the result 
expected in the pure magnetic field scenario. Case III is interesting and quite relevant to real world experiments. 
Compared to the situation in a pure magnetic field, the geometric phase in the presence of a strong tensor Stark 

/ u z B z \ 2 

splitting is significantly suppressed, by a factor y — j which could be ~ 10~ 6 in experiments using polar molecules 

to search for the electron electric dipole moment jTTJ [T^l HH] • 

The expressions listed in Table [TT] are special cases derived from the more general Equation ( 38 1 for the energy 
shift between |1,±1). A priori it is not intuitively evident how Equation (38 1 or the limiting cases in Table P can be 
related to a solid angle in the parameter space of the Hamiltonian. But reformulating the problem in terms of energy 
level shifts allows us to exchange a convoluted geometric construction for a simple algebraic calculation. 



Case Limiting condition Geometric phase difference, A<j} g 



I |Ai|<2B< Alz S z 



uj\T 



ii |a 1 |« m ^,«2b (j^t) 2uj±t 

II n z B z « |Ai| « 2B {^f) ^ T 



TABLE II: Analytical expressions in some limiting conditions, for the geometric phase in a static electric field £ z with a rotating 
transverse magnetic field component B±. In each case, the geometric phase between |1, ±1} up to first order in u± is shown. 
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VI. SUMMARY 



We have worked out some simple examples where the geometric phase arises due to magnetic or electric fields 
whose direction changes over time. Identifying the geometric phase as nothing other than the phase due to off- 
resonant energy shifts enables the use of textbook quantum mechanics to calculate it in every case. In contrast to 
the usual geometric approach, geometric phases due to non-adiabatic evolution and/or non-cyclic paths in parameter 
space can be calculated in a straightforward way using energy shifts. Furthermore, the energy shift approach is 
convenient for understanding geometric phases when complex level structures or multi-dimensional parameter spaces 
are involved. Instead of tracking all the cigenstates as the Hamiltonian evolves along a path in its parameter space, 
the problem is reduced to one of simple and intuitive algebra. This simplification enables the precise calculation of 
geometric phase effects in experimentally relevant situations, where the evolution of the fields experienced by atoms 
or molecules can be quite complicated. Under conditions where the transverse fields are sufficiently small, using 
perturbation theory to calculate energy shifts (as described here) gives particularly simple expressions even when the 
relation to the actual geometry of the situation is unclear. 

The analysis in terms of energy shifts outlined here is not meant to replace the geometric approach. The beautiful 
and general results known from geometric arguments emphasize some of the topological aspects of geometric phases, 
which cannot be easily obtained in the energy shift formulation. However, it is useful to realize that the full machinery 
of the geometric formalism is not in general necessary to calculate the relevant phases. Moreover, from the pedagogical 
point of view, we hope that recasting the geometric phase in terms of familiar concepts (phase differences due to energy 
shifts) may make it seem a little less mysterious. 
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